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The ttN — i> rjN and rjN —* rjN partial wave T-matrices for the eight lowest partial 
waves have been obtained in a three coupled channel model with unitarity mani- 
festly imposed. The two physical channels are vrA^ and rjN, and the third channel, 
ttitN is an effective, but unphysical two body channel which represents all remaining 
processes. The irN elastic phase shifts and the weighted data base of the vrA^ — > r]N 
total and differential cross sections are chosen as the input for the fitting proce- 
dure. A model containing a single resonance in each of the three partial waves that 
dominates the rj production at lower energies is compared with previous analyses, 
based on similar assumptions. A multiresonance coupled channel model is intro- 
duced which significantly improves the agreement with all input data. Our results 
are compared with a complementary multiresonance coupled channel analysis that 
is constrained with elastic and continuum production channels. The inclusion of the 
fourth Pii resonance in the 1440 - 2200 MeV region further improves the agreement 
between the analysis and the data. 



I. INTRODUCTION 



The cross section for //-meson production by pions 

TT^p — s> rin (1) 

has a large peak at an energy close to the rj production threshold, see Figs, la 
and lb. The maximum value, o'toti'^ P V^) = 2.8 mb, is 6 % of o'toti'^ P —* all) 
[|^]. This large peak is usually associated with the Sii(1535) isobar known from vrA^ 
elastic scattering and photoproduction. Recently, Hohler has drawn attention 
to a peculiar feature of this Sn resonance, namely, that the speed plot analysis of 
S-wave ttN scattering has a sharp spike at the opening of the rj channel and no clear, 
independent indication of a resonance in the energy region around 1535 MeV. This 
suggests a strong interplay between the cusp associated with the opening of the rj 
channel and the excitation of the Sn resonance. 

Some time ago Bhalerao and Liu 0] analyzed the then available vrA^ — » r]N data 
in a coupled-channel, single resonance separable interaction model in which the 
reaction proceeds via the formation of A and A^* isobars - and concluded that the 
S-wave T]N interaction is attractive. They found the value of the S-wave scattering 
length of ar,N = (0.27 + zO.22) fm. Arima et al |1 obtained a^TV = (0.98 + iiO.37) 
fm, while Wilkin, based on an S-wave threshold enhancement calculation, quotes 
the value of o^at = (0.55 ± 0.20 + i 0.30) fm 0. The large spread in these values of 
the fundamental QrjN parameter illustrates the need for better understanding of the 
TjN system at low energies. The first step to do this is to obtain a reliable set of 
ttN TjN T-matrices. 

The indication of strong and attractive rjN interaction has led to a speculation 
about the existence of a new type of nuclear matter, namely, quasi-bound ?7-mesic 
nuclei Q . The properties of this new matter are determined by the TjN interaction 
at low energies. 

Good data on rj production in vr^p interaction is missing. The dominant con- 
tribution to the surprisingly big r] production channel is coming from the Sii(1535) 
resonance, the contribution of the Pii(1440) and Di3(1520) resonances is important, 
but not completely clarified. The role of other resonances, even in these partial 
waves, is not at all discussed because of the single resonance character of the model. 
Recently, accurate r] photoproduction data have been obtained by TAPS at MAMI 
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[0] up to = 790 MeV. These data indicate that the D13 resonance contribution 
is smalL 

The objective of this paper is to furnish a set of ttN —>■ rjN partial wave T 
matrices that describe the available data in a straightforward way. It is essential to 
get reliable information on vrA^ rjN on shell T matrices to be able to calculate the 
higher order processes where the 77 production vertex is a part of the higher order 
diagram, for instance in pp pprj. Several publications, dealing with higher order 
processes, have either used the dominant Sn partial wave [p|-[TO[| or have included 



other partial waves 0,|ril] for obtaining the elementary -nN rjN amplitudes. It 



is clear from Fig. la and lb, where a comparison of different model predictions 
with experimental data is given, that only one s-wave resonance is not sufficient, 
the Sii(1535) resonance accounts only for the part of the 77 production total cross 
section, for the energy range of the Sii(1535) dominance. A single resonance model 
for the s-wave only is, therefore, incapable of describing all of 1] production data, but 
it may suffice for the region of the Sii(1535) dominance, which corresponds to p,r less 
than 850 MeV/c. To set the stage, we have tested a simple single resonance model 
for all three dominant partial waves; Sn, Pn and D13 without any background terms 
added. The model describes the large near threshold peak, but it fails miserably in 
comparison with data above T^r = 800 MeV. The recent single resonance model by 
Bennhold and Tanabe (B-T) []T^ is, anyhow, limited to T^r less than 700 MeV. It 
describes the dominant peak pretty well on the gross scale, but significantly fails in 
giving details like the exact peak position, etc. (see Fig. lb). The reason for the 
failure is that the B-T analysis relies on the data of Ref. which suffer from a 
serious beam momentum calibration error We have also used data of Ref. ||13|| , 
but with caution, as elaborated in detail later in the text. The comparison between 
the two models is shown in Fig. la and in the expanded scale in Fig. lb. Considering 
that the model of B-T |1^ used erroneous data of Ref. [0 and has an energy 
dependent form factor, which is not included in our study, the agreement between 
the predictions of our single resonance model and the B-T model, is satisfactory. Of 
course, both single resonance models miserably fail to fit the data above 800 MeV. 

To obtain a better description of the input data we propose a multiresonance cou- 
pled channel model with a smooth background added in a unitary way similar as 
in the Karlsruhe-Helsinki partial wave analysis (K-H PWA) of vrA^ elastic scattering 



16]. A data base consisting of nN K-H PWA of Ref. with the addition of total 
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and differential vrA^ — > r]N cross sections ||1| J13|J17H 23| has been chosen in order to 
perform the fitting procedure of the vriV elastic partial waves up to T.^ = 2727 MeV 
with the additional weighting factors based on the analyses of the data reliability 



by Clajus and Nefkens |T^. Results for the obtained resonance parameters slightly 
deviate from the vrA^ elastic case, but they give acceptable vrA^ — > rjN partial wave 
T-matrices and at the same time predict the T matrices for the rjN elastic pro- 
cess. Manley and Salesky (M-S) |^ have used a multichannel and multiresonance 
K-matrix approach to the coupled channel inelastic ttN scattering with the main 
inelastic channel beeing the continuum production. The rj production process is only 
added in two partial waves, the Sn and F17. The only purpose for including the r] 
production channel is to maintain unitarity when the observed loss of flux can not 
be attributed to any other process. The so called "77 production channel" have to 
account for ALL missing inelasticity within the partial wave and cannot be related 
just to 1] production. Therefore, the prediction for the two r] production partial 
waves T-matrices Sn and F17, which come from the M-S model p^ , p5| should be 
used with care. As the complete study of the rj production process requires all partial 
waves we have made a multiresonance three coupled channel fit with the rjN channel 
explicitly included and used as a constraint. The multiresonance M-S model [p^ , p5 



is complementary to ours in the sense that our results for the resonance branching 
ratios for the tt'^N channel can be compared to the continuum production branching 
ratios in M-S model, and they should (but only roughly) correspond. 



II. FORMULATION OF THE MODEL 

A. The three body coupled channel formalism 

The ttN rjN process is given by the invariant amplitude 

A{w,cose*)+ 4r,B{w,cose*) 

with the standard on-shell partial waves decomposition of A and B: 

in 

qfq;-" { 1=0 



A{W, COS e*)= , <^ yTi+ J(E* +m){E*f + m){W -m)P; (cose*) 



+ J{E* - m){E} - m){W + m)P/+i(cosr; 
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-^Ti-[ ^{E* + m){E} + m){W - m)P/+i (cos r 
1=1 

+ ^{E* - m){E} - m){W + m) Pi {cos 9* 



(2) 



B{W, cos 9*) = 



47r 



oo 
1=0 



E* + m){E} + m)Pl{cos9*) 



-^{Et-m){E}-m)P;^,{cos9*) 

oo 

+ J2Ti-[^{E* + m){E} + m)Pl^,icos9*) 
1=1 

-^{E*-m){E}-m)P;{cos9*) 

W is the total cm. energy, 
9* is the cm. scattering angle, 

q* and q* are the initial pion and final rj cm. momenta, 
E* and EJ: are the initial and final nucleon cm. energies, 
Pl{z) are derivatives of Legendre polynomials, 
Ti± are the nN — > riN T-matrices, 
and m is the nucleon mass. 

The ttN — > rjN Ti+^i- matrices are matrix elements of the three channel partial 
wave T"^^ matrix which is given as: 



/ rjiJL rpJL rpJL \ 
TTTT TTT) nw^ 



T 



JL 



rpJL rpJL rpJL 
rpJL rpJL rpJL 



J 



where various channels are denoted by the index tt for ttA^, 77 for 77 A^" and tt^ for all 
other channels (ttA, pN, nnN, ...). The third channel is effectively described as a 
two body process tt^N with tt^ being a quasiparticle with a different mass chosen 
for each partial wave. We have fixed the channel masses, for each partial wave 
independently. 
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B. Single resonance model 



The simplest possible model to represent the iiN vertex, and which is directly 
comparable to earlier analyses P, p^p!2| , p6| , has been constructed. Only one reso- 



nance without any background terms is used to describe each of the three important 
partial waves Sn, Pn and D13. 

The elastic T-matrix for each resonance is defined as: 

^"'(""' = -w-^lmi2 ' ' = ^- "• 

The partial widths are given by: 

2L+1 

^ ^ for Qc < Qo^ 



\ 2L+1 

and Qc is the cm. momentum of the channel meson c: vr, r], or vr^: 



(4) 



?c = gc(w^) = ^^1^ (5) 

for TTic = "n^Tr, rrirj, m^2. ggj" = qc{M^^) is the cm. momentum of the channel meson 
c at the resonance mass M"^^. 

The exception is Pn because its mass is below rj threshold. For that case we use: 

„ I below riN threshold 

rf"(iy) = i ^ ^3 ' (6) 

^ [100 (^) above r]N threshold 

This definition of partial widths gives the correct threshold behavior for the T- 
matrix. The total widths are given as: 

^ittiw) = Ti\W) + rf (W^) + Ti^{W). (7) 

The inelastic T-matrix is given by 



Tii^m = v/r,f (l^)T,f (W^) c„C2 = n,v,7,', (8) 
and the T-matrix is defined in such a way as to give a unitary S matrix: 

^ = 1 + 2iT, S^S = SS^ = 1. (9) 
5 



Reliable branching ratios of various resonances to the rj channel unfortunately are 
not available except for Sn. As can be seen from Table 1, all values used in our 
calculation are within range of the accepted values [|T|. However, the agreement of 
the elastic vrA^ partial wave T-matrices values with the Karlsruhe - Helsinki phase 
shift analysis (K-H) |T6[ is not entirely satisfactory for the lower energy range T^r < 



700 MeV, and the single resonance model completely misses all higher energy partial 
waves in the elastic channels. The results are shown in Figs. 2a and 2b. The dashed- 
dotted line shows our single resonance model and the dotted line is the B-T single 



resonance model of Ref. |T^ which has an additional energy dependent form factor. 
The B-T model is limited to lower energies. The dots are the elastic K-H vrA^ partial 
wave analyses fl^ . As the K-H PWA does not give the error analysis for the partial 



wave T-matrices in [T^, and the errors are essential to define the statistical weight 
of the analyses, we have identified the errors of the used data in the standard 
analysis as: 

/ W- — thresh \ 

Af^ = 0.005 + I 0.01 + 0.0015 — j {T^^J 

A = 1 GeV 

Wi is the total cm. energy 
thresh ^j^^ total energy at n nucleon threshold 

l^maxl is maximal value of the T — matrix in the chosen energy range. 

The energy range extents up to 2.5 GeV of the total cm. energy. 

The statistical weight in the function is defined in a standard way as well: 

1 



The introduced energy dependence of the statistical weight is inspired by the 



energy dependence of the error analysis of Ref. [|T5|. It steadily raises with energy, 
but does not exceed the value of 0.02 in the units of Ref. ]16| . 

The total cross sections for the rj production are compared with the single reso- 
nance model (dash-dotted line) of this article, and the B-T single resonance model 
of Ref. [|l2| (dotted line) in Figs, la and lb. The experimental data for the rj pro- 
duction differential cross sections are compared with both models in Figs. 3a-3d. 
The full and dashed lines in Figs. 1-3 will be explained later. 
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The failure of single resonance model to describe the experimental facts is not 
unexpected. Even at the lower energies, background terms which are not included in 
the single resonance model are of some importance, and their omission can account 
for the observed discrepancy with the input. The comparison of the parameter 
free single resonance model, of this article (dash-dotted line), with the similar B-T 
model (dotted line) shows a fair likeness in the energy range where the later 
analysis has been given. The differences between the two models are attributed to 
the different input data base and to the additional energy dependent form factor used 
in [O. However, as is seen in Figs. 1-3, both models show reasonable agreement 



with the input data for the important nN elastic partial waves Sn, Pn and D13 at 
lower energies. Our analysis predicts strong deviations from the fitted vrA^ elastic 
partial waves at higher energies, for the Sn partial wave in particular, so it is to 
be expected that the results of [|1^ suffer from a similar disease for higher energies. 
Similar statements can be and are valid for the total and differential cross sections 
of 7] production (Figs. 1 and 3). The phenomenon observed is consistent with the 
fact that the single resonance model covers only the energy range of the first peak 
in ati^Tip T]n). 

Let us mention that in the described single resonance formalism it is not possible 
to extrapolate the model to include more than one resonance per partial wave in a 
straightforward manner without directly violating the S-matrix unitarity. 



C. A unitary multiresonance model 

In order to fit the ttN elastic amplitudes of Ref. [|l^ better then in single reso- 



nance model we have introduced a manifestly unitary model that enables including 
more than one resonance and background term per partial wave. It is constructed 
following the commonly accepted method developed in []T3|, and originally used in 
TiN elastic partial wave analysis (PWA) |]rB|. ^ 

As the 1] meson is a pseudoscalar, isospin zero particle, it mixes by isospin vio- 
lation with the TT^. We have chosen the following three coupled channels to set up 



*The simple recipe for the modification of the S-matrix, which consists in directly adding 
extra resonances and a smooth background is manifestly nonunitary. 
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the model: the vrA^, rjN and a third, an effective two body channel labeled vr^A^, 
which inclusively contains and represents all remaining, even three body channels 
(ttA, pN , ttttN, etc.). The objective of the procedure is to simultaneously achieve 
a good representation of the input vrA^ elastic T-matrices, and the experimental 77 
production data (total and differential cross sections) by the values coming out of 
the model. 

We distinguish three basic steps of the presented analysis. 
Step 1 - formalism 

is given as: 



The multichannel T matrix taken over from 15 



^af = E //^(^)v^7„tG?(5)7?V^//^(5) (10) 

a, 6 = TT, rj, TT^. The initial and final channels couple through intermediate "particles" 
or resonances labeled with i and j, and 

fa\s) = I '-^—=\ (11) 



Pais) = ^ (12) 

where s = and is the meson momentum for any of the three channels given 
as: 

^ - (m + maY) (W^ - (m - m,)^) 
Qa = qa{W) = ^ — . (13) 

The mass parameter 171^^2 for the tt'^N channel is fixed prior to minimization to the 
mass value at which the partial wave inelasticities show the opening of the first 
inelastic channel: 



Partial wave 


Sii 


Pii 


Pl3 


Di3 


Di5 


Fl5 


F17 Gi7 


m^2 (MeV) 


450 


380 


370 


380 


400 


370 


650 450 



7j^^ are free parameters and will be determined by the fitting procedure. For the 
Qia and Q2a parameters we choose the values 

Qitt = Qi-K = fn-iT 

Qlrj = Q2ri = niri (14) 

Qi,2 = Q2^2 = 400 MeV 
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Gf^ is a dressed propagator for partial wave JL and particles i and j; 

Gi^s) = G^\s) + E Gf,\s)i:i^{s)Gil^{s). (15) 

A;,Z=1 



The bare propagator 



(.) = ^ (16) 

^ Si- s 



has a pole at the real value Sj. The sign Cj = ±1 must be chosen to be positive for 
poles above the elastic threshold which correspond to resonances. The nonresonant 
background is described by a function that consists of two terms of the form ([TBI) 
with pole positions below vrA^ threshold. For that case signs of terms are opposite. 
The positive sign correspond to the repulsive and negative to the attractive potential. 
E^/" is the self-energy term for the particle propagator: 

ni{s) = T.^i!:<\sW,j^ (17) 

a 

The imaginary part of (^^J"{s) is the effective phase-space factor for the channel a: 



JL/ 



JL, 



tfM^i^r ,_, ^y) _ (19) 



The real part of ^■^^{s) is calculated using a subtracted dispersion relation 

-So /-oo F^^{s') 

TT Jsa (s'-s)(s'-So 

where = (m + rria)'^. 

The advantage of this approach is that it manifestly maintains the 5* matrix 
unitarity for any number of resonance and/or background terms. The disadvantage 
is that the connection of the parameters ■yfj" and Si with the conventional resonance 
parameters Mf^ and T^J^" is not direct [|T^, but has to be calculated. 
Step 2 - data base and fitting procedure 

The input parameters for the fitting procedure are Sj and which determine the 
bare propagator and self energy term for the particle propagator, see Eqs. (|T^),(|T^) 



and (|TBp respectively. The parameters Qi^ and which occur in the form factor 
given in Eq. ([1^) have been fixed to the mass of the channel a meson. The ones 
subtracted dispersion relation given in Eq. (|T9|) is solved numerically with the 
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subtraction constant sq = Sa and $;[^(so) = 0. The stability of the solution has been 
tested by calculating and reproducing the initial imaginary part, see Eq.(|^). The 
numerical integration has been performed using the adapted Gaussian quadrature 
method with no significant dependence on the number of points. We should mention 
that the dispersion relation has been calculated only once, and tabulated for further 
use to save the CPU because the parameters which form the integrand are not varied 
in the minimization procedure. 
Data base 

As the input data to the minimization procedure we have used: 



1. the K-H partial wave ttN analyses [IG 



2. total cross sections for the vrA^ r]N process |l3| , p!7| -p^,pS 



3. differential cross sections for the vrA^ — > rjN process [ p!^ , p!^^ - |20|p^ , |23|p8| 



We have fitted the ttN elastic T matrices for eight / = 1/2 partial waves: 5*11, 



Pii, Pi3, Di3, Di5, Fi5, Fn and Gn using the phase shift analysis K-H [16 at 



90 energies from threshold to 2.5 GeV of the total cm. energy. We have also 



tested the use of (CMU-LBL) nN elastic partial wave analyses |15] without any 
notable differences. At the inception of this work the only PWA that covered the 
high energy was K-H [|16|. Since then the VPI group [27| has extended its analysis 



approximately to 2.1 GeV. We do not anticipate that the use of VPI PWA would 
introduce any substantial changes into the conclusions coming out of the present 
analysis. Anyhow, a new analysis based on the VPI PWA is planned to be the 
subject of a future research. Finally, let us comment that we should have used 
the total set of experimental data for vrA^ elastic processes instead of limiting the 
input to the model dependent information coming from PWAs. However, we have 
assumed that existing PWAs represent the data adequately enough so we can avoid 
an enormous CPU time consumption. 

The data for the second, riN channel, are the vrA^ — > rjN differential cross sections 
at 81 energies and total rjN production cross sections at 67 energies |I3| , p!7| - pB|J28| . 
The statistical weights of for the rjN data, used in the minimization procedure, have 
been based on the analysis of the world data give in Ref. |]13. In some cases the 



published statistical weight had to be modified. Problems of consistency among 
different measurements have been extensively discussed in WM, and the discrepancy 
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for the lower energies of t] production differential cross section of data of Ref . [|1^ has 



been claimed. It has been argued in Ref. the data of Ref. ||13| are systematically 
too low, and that it is due to an error in the beam momentum calibration which 
makes the data at lower energies essentially unusable even if one tries to correct them 
by a momentum shift. The direct reason for that is a strong momentum dependence 
of the lab -vv- cm. transformation Jacobian. However, at higher energies, it is safe 
just to perform a 4 % momentum shift downward. Therefore, the systematic error, 
in addition to the published, statistical one, has been added to the questionable data 
sets reducing their statistical impact. The statistical weight of all differential cross 
sections of Ref. [|13| has been reduced by the factor of 2. In addition, the statistical 
weight at lowest energies oi W = 1511 MeV, 1542 MeV and 1571 MeV has been 
even more reduced so that their importance in the analysis is practically eliminated. 

A systematic error of 0.01 mb/sr has been added to all rj production differential 
cross sections. We have decided to do so because of the fact that the quoted errors 
in all measurements have been of statistical origin only, so most of the errors have 
been unrealistically low. 

The data taken at the energy of 1507 MeV in Ref. [|l^ as well tend to be too low 
when the total cross section is calculated [Q and compared to the "world trend". 
Therefore, the additional systematic error of 0.04 mb/sr have been added to these 
data. 

For similar general reasons the systematic error of 0.02 mb/sr has been added to 
all the data of Ref. [EH . 

All "two star" cross sections, as given in Ref. are taken with the increased 
statistical factor of 10, while the statistical factor of the remaining cross section data 
is kept to be 1. 

The direct consequence of our choice of the data base statistical weights is seen 
in Fig. 3a: Ref. reproduces the data set of Ref. [|13| at lower energies ( W 
= 1511, 1546 and 1571 MeV) while our result tends to reproduce other sets much 
better. Additional, precise measurements for that process are needed in order to 
eliminate the present uncertainties. 
Fitting procedure 

We have fitted the 8 above listed lowest 1 = 1/2 partial waves in the following 
manner: 
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The number of resonances and the shape of the background is determined by the 
choice of bare propagator and self energy term parameters. We have decided 
to start from the number of resonances in each partial wave as given in P|JIB|]; 
the background is represented by two resonant functions with the constraint 
that the pole position is far outside the physical region. Only one resonance in 
the Pi3 and F15 partial waves has been used in our analysis, higher (uncertain) 
resonances have been dropped. The existence of the second P13 resonance has 
been recently suggested, and it was used in the M-S t^N analyses p4| , p5 . 



b. In general, we need up to two background terms. They have been numerically 

represented as tails of resonances having their pole positions far outside the 
analyzed energy rangeQ. 

c. Error analysis for resonance parameters has been done on the basis of MINUIT, 

imposing the confidence level of 70 9l 



We have used a standard MINUIT program using as much as 132 parameters in the 
final run. 

The minimization has been complicated, and a lot of technical tricks had to be 
used to avoid occurrence of local minima which are hard to handle for a minimization 
with such a number of parameters. 

The result of the fitting procedure gives us a full three channel T-matrix, with 
submatrices describing the vrA^ — > rjN and rjN rjN processes. 
Step 3 - resonance parameters 

The pole positions, resonance masses and widths have to be obtained numerically 
from the full partial wave T-matrix, following directly the technique developed in 



1^ . For the convenience of the reader we shall briefly reproduce the essential steps, 
angular momentum indexes are suppressed. 

The poles of T-matrix given in Eq. (110) are found solving the following equation: 



detG"^ = 0. (20) 
The eigenvector of the matrix H = G^^ has been found at the pole position Spoie: 



*More than one background term was needed in order to obtain some nontrivial, but 
smooth energy dependent behavior of the background terms. 
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Y,H,,MXj = 0. (21) 
j 

We have defined quantities: 

Vc = J2 ^i<^Xi (22) 

i 

which define the coupling of resonance i to the channel c. We consider the width to 
be an energy-dependent quantity involving the phase-space factors: 

C 

Vc = \vf. (23) 
Near the resonance we parametrize the T-matrix as: 

Tat = {Bat - Sat)/2t + ^ F^/%D-%Fy' bH;' (24) 

cd 

where Bat is a background S-matrix and the generahzed Breit-Wigner denominator 
is: 

D{s) = r-s-cY,y,<!>,is). (25) 

c 

The real constants r and c are chosen so that -D(spoie) = 0. The resonance mass, 
width and the branching ratios are defined as: 

Re D{M^) = 

^ Im D{M') 

MRe D'{M^) ^ ' 

_ y.F,{M')T _ 

where D' is the derivative of the generalized Breit-Wigner denominator (pSj). The 
obtained resonance parameter values are given in Tables 2 and 3. 

III. RESULTS OF THE UNITARY MULTIRESONANCE MODEL WITH 

THREE Pii RESONANCES 

We have obtained the partial wave T-matrices for vrA^ elastic scattering (see 
Figs. 2a and 2b), vrA^ — rjN (see Figs. 4a and 4b) and rjN elastic scattering (see 
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Figs. 5a and 5b), on the basis of our fit to the ttN elastic and nN —>■ rjN data. 
The agreement with the input vriV elastic K-H PWA T-matrices is given in Fig. 2, 
while the agreement with the input tcN rjN differential cross sections is given in 
Fig. 3. The full, dash-dotted and dotted lines systematically denote the results of 
our multiresonance, and the two single resonance models, this publication and B-T 
respectively. The T-matrix for rjN elastic scattering is a prediction. All observables, 
obtainable on the basis of our results for the partial wave T-matrix of the vrA^ rjN 
and rjN elastic reactions, with the exception of the nN — > rjN differential cross 
sections, which are input data, are henceforth a prediction. Our T matrices can be 
used as input to rj production calculations in reactions such as pp pprj. 

As can be seen, single resonance models give a reasonable agreement only in the 
region of the Sii(1535). They can be used for higher order calculations only, for a 
limited energy range where one resonance per partial wave dominates. However, the 
analyses presented in this article, which is based on more than one resonance per 
partial wave, can be used over the full energy range. 



The results of the multiresonance model for the ttN interaction by M-S pj] can 
be used as a consistency check. The parameters of both models, for the same number 
of resonances per partial wave ( those given by the PDG [IH) are listed in Table 2. 
We may call these analyses complementary because the inelastic part of the ttN 
partial waves is constrained by two complementary processes: rj production in our 
case, and continuum production in the case of Ref. Therefore, the vr^ part of 
our analysis should roughly correspond to the vrvr part of another analyses, and the 
parameters of the 1] production partial waves explicitly included in other analyses 
should correspond to our findings. Of course, masses and widths of resonances should 
correspond to the values given by the K-H and CMU-LBL nN elastic analyses, which 
are generally accepted by PDG ffj]. 

As can be seen from our Figs. 1-3, the multiresonance model, based on the 
standard number of resonances describes the input vrA^ — > vrA^ and vrA^ rjN data 
fairly well. Of course, the 1] production cross sections are as well correctly described 
in the full energy range. The structure of the irN elastic partial waves is not 
entirely reproduced. The tendency of smoothing elastic partial waves, as it has 
been already indicated previously, exists when the inelastic channels are explicitly 
included. Therefore, we are tempted to conclude that inclusion of inelastic processes 
imposes some restrictions on the elastic channel forcing partial waves to have less 
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structure then in |jT6 . 

Comparing the available analyses we conclude: 



1. The masses and vrA^ elastic branching ratios generally agree for all three PWA, 
with the exception of the Pn. This problem is discussed later. 

2. All partial waves except the Pn show reasonable agreement in both multichan- 
nel models. 

3. The two higher Pn partial wave resonances of M-S and our analyses do not 
agree. The branching ratio for rj production in our model is about 90 %. 
The M-S model predicts almost a 60 % branching ratio to the vrvr chan- 
nel, leaving no freedom for any flux going to the t] production channel. The 
disagreement is obvious, and we shall offer a natural explanation. 



IV. RESULTS OF THE UNITARY MULTIRESONANCE MODEL WITH 

FOUR Pii RESONANCES 



Inspection of the resonance parameters of Ref. reveals the following (see 
Table 2.): 



the total width of the Pii (1440/135) is different from the ones of K-H |jT6[, CMU- 
LBL [0 and Arndt et al. |2j 



the total width of the Pii (1710/120) 



is different from the ones of K-H ||T6| , CMU 



LBL [|T|] and Arndt et al. |2 



• the mass the Pii (2100/200) much lower than the ones of K-H |T6[ and CMU- 
LBL [0 

• the mass of the 013(^330/265) is shifted from 2080 to 1804 MeV 

Therefore, we suspect that a part of the physics in the vicinity of the 1800 MeV 
mass region is not entirely taken into account. 

We assume that there is another degree of freedom in the Pu partial wave, in 
the form of another resonance. 
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This possibility of having four instead of three resonances in the Pn channel leads 
to a fit which is shown in the figures with dashed lines. The resonance parameters 
are given in Table 3. The T-matrix pole positions for K-H [|T^], CMU-LBL | |15| |, 
M-S and our analyses (3 and 4 poles) are given in Fig. 6Q. As can be seen in 
Fig. 6, the pole positions of the three resonances are well established in the classical 
ttN elastic PWAs and are fairly close. The Pn T-matrix pole positions for the 
M-S analysis (given by filled circles) are quite different from the recommended 
values. The three pole version of the Pn of our calculation is also quite compatible 
with the K-H and CMU-LBL pole positions. The four pole version of the Pn shows 
agreement for the lowest pole, the next two poles are near in masses, but distinctly 
separated in the complex energy plane, while the fourth pole is somewhat lower than 
K-H. In our analysis the fourth resonance is strongly inelastic, going mainly to the 
'q production channel, therefore it is not to be expected that either of the K-H or 
M-S analyses could determine it with great precision. 



V. CONCLUSIONS 

1. The addition of another resonance in the Pn partial wave definitely improves 
our fit to elastic and inelastic data in all channels. Various quark models also 
predict four and even five Pn resonances in the energy region of 1440 to 2200 
MeV |30|J31|]. 

2. The changes in vrA^ elastic partial waves are negligible because all the reso- 
nances (with the exception of the first one) are strongly inelastic. 

3. The changes in r] production and rjN elastic channel T-matrices are clearly 
visible in all partial waves. As it is to be expected, the Pn is significantly 
changed, while other partial waves do show some variation. However, let us 
draw the readers attention to the fact that Sn, Pn and D13 partial waves are 
the dominant ones, while the contribution of other channels is at least the 
order of magnitude lower. So, even if the relative change in Figs. 4 and 5. is 



*CMU-LBL and our analyses give the exact T-matrix pole positions, while the poles for 
the remaining analyses have been approximated by M — iV /2. 
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large for other partial waves, the change at the absolute scale is comparably 
small. 



4. The inconsistency problem between the two inelastic PWAs goes away. The 
two Pii resonances in the vicinity of 1750 MeV are responsible for the con- 
tinuum production, this is to be compared with 1717 and 1885 resonances of 
Ref. P^, the third resonance at 2215 almost completely couples to the rjN 



channel, with a very small branching ratio to the continuum production. 

The M-S PWA could not easily see the additional fourth Pn resonance 
as it mostly couples to the r) production channel, and it is not, in their case, 
explicitly included. 



6. The elastic vrA^ analyses 1T6|JT^ have problems of determining the number and 
parameters of the resonances going to inelastic channels. Some changes in the 
resonances with a small elastic branching ratio should be easy to obtain in the 
single channel formalism. 

7. If the Brown data for T^r > 900 MeV are shifted to lower energies by 30 
MeV and data below 900 MeV are omitted then the resonance parameters for 
our solutions containing 3 and 4 resonance in the Pn are shown in Table 4. 
The results are virtually identical (within errors) with the solutions using all 
originally published Brown data . That indicates that our PWA is showing 
robustness with respect to possible errors in input. 

Of course, adding extra Pn resonance is only one way of making our analysis compat- 
ible to M-S PWA, there might be other, equally good explanations of the apparent 
disagreements when the inelastic branching ratios are compared. One technical de- 
tail has to be kept in mind: the background parts are represented by two resonances 
with poles kept far outside the range of interest. Therefore, the background is also 
fitted. The proper way should be to calculate the background in some model (for 



instance a cloudy bag model ), to fix it, and then fit just the resonant part ||32|. The 
interference of background and resonant parts might shift the resonance parameters 
notably. 

The T matrices obtained in this study are an essential ingredient for calculating 
7] production reactions D(7r,NN)?7 and D(7r,NN)7r. The off mass shell extrapolation 
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procedure for the rj production amplitudes can not be determined in this formahsm. 
We do hope to learn something about these effects from the higher order processes, 
assuming that the off-mass shell behavior of the rj particle is similar as of the pions 



*This work has been partly supported by the EC contract CIl*-CT-91-0894 and DOE 
contract DE - FG03 - 88ER40420/A606. 
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TABLE CAPTIONS 



Table 1. Resonance parameters of the single resonance model. 



Table 2. Resonance parameters of the multiresonance model with 3 Pu resonances. 
The results of elastic vrA^ analyses [jl],|l6|,|l5| are given in the first column. The 
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results of the partial wave analysis of this article, as well as results of the PWA 
of Ref. 1^ are given in columns 2-10. 



Table 3. Resonance parameters of the multiresonance model with 4 Pii resonances. 

Table 4. Resonance parameters of the multiresonance models with 3 Pu and 4 Pn 
resonances using the data of Ref. [|13|, for T^- > 900 MeV shifted downward by 
30 MeV. Brown data below 900 MeV are omitted from fitting procedure. 



FIGURE CAPTIONS 

Fig. 1. The total rj production cross sections. The experimental data are taken 
from the literature, as it has been indicated in the figure itself. The dotted 
line is the approximation that only one resonance has a strong branching ratio 
to the rj production channel, the results are taken from [^. The dash-dotted 



line is the result of this calculation, and is based on similar assumptions. 
The full line is the result of the three coupled channel multiresonance model 
presented in this publication with the number of resonances given by the PDG 
m, namely three in the Pn partial wave. The dashed line is the four Pn 
resonance model. The full lines and the dashed lines show notable difference 
for the energies above T^^ > 800 MeV. The four resonance model for Pn is 
closer to experimental data. 

Fig. la is the entire energy range available. 

Fig. lb is the same data, with an expanded scale for T^r up to 800 MeV. 
Fig. 2. The vrA^ elastic partial waves. The filled circles are the result of the single 



channel ttN elastic PWA given in Ref. |jT6[. The used PWA does not give 
the error analyses for the partial wave T- matrices in |T^, so the error bars 
given in the figure are defined in the text and reflect the statistical weight 
of the data set used in the minimization procedure. The dotted line is the 
approximation that only one resonance has a strong branching ratio to the 



rj production channel, the results are from |12l|. The dash-dotted line is the 
result of this calculation, and is based on similar assumptions. The full line is 
the result of the three coupled channel multiresonance model presented here 
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with the number of resonances given by the PDG 1^, namely three in Pn 
partial wave. The dashed line is the four Pn resonance model. The full and 
the dashed lines are practically indistinguishable for all elastic partial waves. 

Fig. 3. The comparison of the t] production differential cross sections with the 
afore mentioned PWAs. The experimental data taken from the literature are 
defined in the figure. Curves are as in Fig. 2. 

Fig. 4. The partial wave T matrix for the t] production channel. The dash-dotted 
line is the single resonance approximation for Sn, Pn and D13 partial waves 
of this article. The full line is the result of the three coupled channel mul- 
tiresonance model presented here, with the number of resonances given by the 
PDG |l|], namely three in Pn partial wave. The long dashed line represents 
the increase to four resonances in the Pn partial wave. Differences between 
long dashed and full lines are quite notable in this figure. 

Fig. 5. The partial wave T matrix for the riN elastic channel. Differences be- 
tween dashed (4 Pn resonances) and full lines (3 Pn resonances) are quite 
pronounced. 

Fig. 6. The pole positions for all four partial wave analyses, mentioned in this 
article. The notation is given explicitly in the figure. 
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Table 1 



Resonance P.D.G parameters Used parameters 

811(1535) m = 1520 to 1555 1535 

r = 100 to 250 150 

= 0.35 to 0.55 0.40 

= 0.30 to 0.50 0.42 

Pll(1440) m = 1430 to 1470 1440 

r = 250 to 450 250 

= 0.60 to 0.70 0.60 

= 0*) 

D13(1520) m = 1515 to 1530 1520 

r = 110 to 135 120 

x^ = 0.50 to 0.60 0.55 

Xr, ~ 0.001 0.0015 

*^Xjj — i.e. = is used below, = 100(^)^ above rjN threshold (see text). 
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